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Motivated by recent experiments on biomimetic membranes exposed to several aqueous phases, we theoretically study the mor- 
phology of a membrane in contact with a liquid droplet formed via aqueous phase separation. We concentrate on membranes with 
negligible spontaneous curvature. At small droplet volume, bending energy dominates and the droplet is only partially wrapped 
by the membrane. At large volume, this configuration can become unstable and undergo a discontinuous transition to a state, 
in which the droplet is (almost) completely wrapped by the membrane. A morphology diagram, showing the parameter region 
where such budding transition occurs, is constructed as a function of the membrane tension and the intrinsic contact angle of the 
liquid with the membrane. The effects of spontaneous curvature are discussed qualitatively. 
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1 Introduction 

Recently biomimetic membranes exposed to several aqueous 
phases have been introduced experimentally^^. They are 
found to exhibit a number of interesting and surprising phe- 
nomena, such as partial to complete wetting transitions ^ , bud- 
ding'^, and membrane tube formation^, which are not yet fully 
understood. The origin of these diverse phenomena is the 
competition between the bending rigidity of the membrane 
and the interfacial tensions of the participating phases. Es- 
sentially we are dealing with wetting phenomena on surfaces 
I which are flexible and can attain many different morphologies. 

A particularly interesting process we shall focus on here 
is droplet-induced budding which represents a morphologi- 
cal transition from a state where the liquid droplet is partially 
wrapped by the membrane to a state where the droplet is al- 
most completely wrapped by the membrane. The resulting 
bud is connected to the original membrane by a small neck, as 
shown in Fig. [T] 

Budding is an important and frequent cellular process. For 
example, it represents an important step during endo- and ex- 
ocytosis of all membranes. Endocytosis leads to the formation 
of transport vesicles^, which allow communication and trans- 
port of biomolecules between different organelles. Likewise, 
budding also occurs during viral replication processes^. In 
addition, there have been proposals for using synthetic mem- 
branes for technological applications, e.g.^ ^. In such a case, 
it is often desirable to mimic real cellular processes, includ- 
ing budding, as the modus operandi. This is a particularly 
promising approach in the context of microfluidics, the minia- 
turization of fluidic operations. 
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Fig. 1 Droplet-induced budding transition: The membrane separates 
the exterior aqueous phase 7 from the two coexisting interior 
aqueous phases a and p . The a droplet partially wets the 
membrane. 



Given the relevance of budding for biological and 
biomimetic systems, it is important to understand its possi- 
ble mechanisms. Budding can be induced by intramembrane 
domains as first predicted theoreticalhM^^ and confirmed ex- 
perimentally by optical microscopy^^^^. The budding pro- 
cess then depends on the elastic properties of the membrane 
domains and on the line tension of the domain boundary. For 
biological membranes, the domains may contain assemblies 
of proteins, which are inserted into the membrane with a 
prefered orientation and, thus, induce a spontaneous curvature 
of the domain. 

In this article, we consider a novel budding mechanism 
where the driving force is the interfacial tension between the 
aqueous phases, see Fig. [T] Droplet-induced budding is some- 
what similar to the encapsulation of solid or rigid particles, 
to which the membrane adheresB^ElJ^ However, in contrast 
to such particles, the liquid droplets considered here change 
their shape during the budding process. While such a bud- 
ding phenomenon has been reported experimentally^, there is 
no detailed theoretical analysis yet. It is our aim here to shed 
light on this new budding mechanism. In particular, we ad- 
dress the required recipes for budding to occur and whether 
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Fig. 2 (a) Geometry of an axisymmetric a droplet. The coordinate 
along the symmetry axis is denoted by z, the distance from this axis 
by p, the arc length by s, and the local tilt angle by All these 
variables depend on the contour parameter t. The contact line is 
located at ^ = . The gray area corresponds to a spherical cap with 
tangent angle 0. (b) Enlarged view close to the contact line: 
Intrinsic contact angle between the two planes that are tangential 
to the aj8 interface and to the smoothly curved vesicle membrane, 
respectively. 



the transition is continuous or discontinuous. We shall focus 
on the limiting, yet instructive, case, for which the bud size is 
small compared to the original membrane area and the mem- 
brane spontaneous curvature is negligible. Our calculations 
show that for a sufficiently low membrane tension, there is 
a critical droplet volume beyond which the partially wrapped 
configuration is unstable and a bud is formed. How low the 
membrane tension must be for a budding transition to occur 
depends on the wetting properties of the liquid on the mem- 
brane. This dependence can be summarized in a morphology 
diagram. Finally, we discuss how the membrane spontaneous 
curvature may modify the budding transition. 



2 Theoretical Description. 

The starting point of our calculation is the theory described 
for membranes and vesicles in contact with two aque- 
ous phases. Here, we consider the geometry shown in Fig. [2j 
The total energy E of the latter system has several contribu- 
tions. First, it contains a term that depends on the volume of 
the liquid droplet and its pressure difference Pya = — 
with the exterior phase 7. The aqueous phase j3 is taken to be 
much larger than the droplet a, which implies Fy — P^ = 0, 
and the j3 phase volume contribution can be ignored. Second, 
the aj3 interface with area A^^jg and interfacial tension I^af5^ 
well as the two membrane segments with areas Aaj and Aj^y 
and membrane tensions I^ay and L^gy contribute three surface 
terms. The membrane tensions T^ay and Z^g^ are mechanical 
tensions that can be related to the anisotropic pressure tensor 
across the membrane. Thus, the tension Z^g^ is well-defined 
and remains, in general, nonzero in the limit in which the two 
aqueous phases j8 and 7 become identical, seeP. In princi- 
ple, the contact line provides an energetic contribution pro- 



portional to its length L(^pY and to its line tension A, but for 
simplicity we have neglected this contribution. Finally, the 
vesicle energy contains the bending energies E'be,; of the two 
membrane segments, with / = a,j3. The bending rigidity, the 
mean curvature, and the spontaneous curvature of the mem- 
brane are denoted by K, M and m respectively. We further 
assume that the bending rigidities and sponteneous curvatures 
are the same in both membrane segments. The Gaussian cur- 
vature term is neglected since we will not consider any topo- 
logical change. Thus, the total energy has the form 



E= ^ JdAiy[2K{M-mf^i:iy]^ZapAo,p^PyaVa- 



(1) 



As shown in Ref .^2, the force balance along the contact line 
is characterized by the intrinsic contact angle Oi^ of the liquids 
at the membrane which satisfies 



-a/3 



COsOin, 



(2) 



and represents a hidden material property of the system. The 
definition of the intrinsic contact angle is shown schematically 
in Fig. [2|b). It represents the contact angle at the nanome- 
ter scale between the aj3 interface and the membrane sur- 
face. Substituting Eq. ^ into ([T]) and defining the rescaled 
areas Aiy = A/y(Zc^^/K-) and Aaf^ = A^p{i:^i^/K), volume 
Va = Va / K*) , membrane curvature M = M{k /ILfy^^ ) ^ , 
spontaneous curvature fh = m(/c/Eofjg)^/^, pressure 



ya 



Pya ( ^/^lcl5 ) membrane tension Z^y = Ziy/Zo^j^ , we ob- 

tain the rescaled energy 

E = E/k = J dAiy2{M-jn)^ ^ cos GinAay^t^pyA^Q 



^af5 +^7a^a7 



(3) 



where Ame = A^y + A^g^ is the rescaled total membrane area. 

We consider axisynmi etric shapes which we calculate using 
the same procedure as in^^ We choose the symmetry axis 
to be the z-axis. The distance from this axis will be denoted 
by p . The vesicle shape is then uniquely described by its one- 
dimensional contour {z{t)^p{t)), as shown in Fig. [2ja). The 
parameter t is the contour parameter, which varies over two 
fixed intervals: to = < t < ti and ti <t<t2, corresponding 
to the two membrane segments Aay and A^g^. Here ti is the 
position of the contact line, and we will focus on the limit of 
large t2, and thus large Aj^y. Two additional quantities that play 
an important role in the theory are the arc length s{t) and the 
tilt angle V^(^). Using this parameterization, the total energy E 
of the vesicle can be written as 



E 

271 
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where the tangent angle 6 corresponds to the angle between 
the ajS interface and the horizontal plane, see Fig. [2ja), the 
rescaled curvature radius R^^^ is determined by the Laplace 
equation Pya = and the Lagrange functions 



\ , I w' sm i/A ^ 1 - , 



s' p 

2 J 



Pyi p^s' sin 1/A + T(p' - / cos . 



(5) 



The Lagrange multiplier T is used to ensure the geometrical 
relation p' = s^cosij/. The primes correspond to derivatives 
with respect to contour parameter t. 

The first variation of the energy E along the membrane sur- 
face leads to the following Euler- Lag range or shape equations 



cos y/smy/ 



:COSl/A+ ^pcosl/A+ ^sin Xj/, 
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sin^ y/ 
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(6) 



while its first variation along the contact line leads to the 
boundary conditions 



(7) 



with the intrinsic contact angle 0in as given by Eq. The 
overdots now denote derivatives with respect to the arc length 
s, rather than the contour parameter t. Furthermore, we impose 
the condition that the membrane is essentially flat for large t2. 
In practice, this is done by imposing = \ = \ \(r\ < £ at 
finite, but very large ^2- We typically choose £ to be of order 
10~^ and t2 to be of order one million of the discretization 
steps. These differential equations are then solved using the 
standard fourth order Runge-Kutta method—. 

3 The budding transition. 

In this paper, we concentrate mainly on membranes with neg- 
ligible sponatenous curvature, m = 0. In this case, the mor- 
phology of the system is determined by three independent di- 
mensionless parameters: the intrinsic contact angle 0in, the 
membrane tension JLj^y and the droplet volume V^. We note 
that generally the membrane tensions JLay ^ £/3y. They are re- 
lated by the equation for the intrinsic contact angle, Eq. ([2]). 
The effects of non-zero spontaneous curvature will be briefly 
discussed in section 6. 

Fig. |3ja) shows the typical energy curve as a function of 
volume for a given intrinsic contact angle and (low) membrane 
tension. As a representative example, we have taken = 45°. 
In the calculations, we have also used a flat membrane with 
no liquid droplet as the reference surface. As mentioned, we 



consider the limit in which the area of this reference surface 
becomes large. Thus the energies shown in Fig. [3] are the de- 
viations from the energy of this flat reference surface. Initially 
when the liquid volume is small, bending energies dominate 
and the membrane bends very weakly, see Fig. [Sj^e). The 
system thus behaves in a similar fashion to the usual wetting 
geometry. With increasing liquid volume, the interfacial en- 
ergies become more important and compete with the bending 
terms. A rather useful naive concept to have in mind is that 
the interfacial energy terms scale as the square of the length 
scale, while the bending term has no explicit length scale de- 
pendence. At higher volume, as shown in Fig. [3jf), while 
the membrane is highly deformed, the system is still in the 
partially wrapped configuration. This remains the case until 
a certain critical volume above which it becomes unsta- 
ble and assumes the completely wrapped configuration (Fig. 
[3jg)). From here on, increasing the liquid volume further in- 
creases the size of the spherical bud and at the same time re- 
duces the size of the neck (Fig. |3jh)). 

It is important to realize that the partially and completely 
wrapped configurations may coexist over a range of volumes. 
As a result, the budding transition described above exhibits 
a hysteretic behaviour. The threshold volume at which the 
morphological instability occurs depends on whether we are 
increasing (V^) or decreasing (V2) the liquid volume. Be- 
tween and V2, there is an energy barrier for a morpho- 
logical transition between partially and completely wrapped 
configurations. In this volume range, the shape equations ^ 
and ^ have three solutions. The two lower energy configura- 
tions correspond to partially and completely wrapped config- 
urations, while the third provides the energy barrier. The latter 
is plotted as dotted lines in Fig. [3] 

Increasing the membrane tension. Fig. [3jb), we find that it 
becomes harder to bend the membrane. Thus, the budding 
transition occurs at a higher droplet volume. Furthermore, 
the volume range, for which both the partially and completely 
wrapped configurations coexist, increases with increasing ten- 
sion. In fact, if we increase the membrane tension even further. 
Figs. [3jc-d), we find that the two configurations effectively 
always coexist. In Fig. [3jc), the completely wrapped configu- 
ration eventually becomes the global minimum configuration 
for large volumes. In contrast, the partially wrapped configu- 
ration represents the state of lowest energy for all volumes in 
Fig.lld). 

4 Morphology diagram. 

The behaviour described in the previous section is, in fact, 
typical for all values of the intrinsic contact angle. Thus, we 
may construct a morphology diagram as a function of intrin- 
sic contact angle and membrane tension. As shown in Fig. 
|4j there are two regimes: (I) a low tension regime, in which 
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Fig. 3 Total energy of the system as a function of the droplet 
volume for various membrane tensions. At low tensions (a-b), the 
partially wrapped droplet (full lines) become unstable when the 
droplet volume V exceeds the threshold volume , whereas the 
completely wrapped droplet (dashed lines) becomes unstable for 
V < ^2*- The dotted lines correspond to energy barriers. In (c), the 
completely wrapped configuration is only slightly more preferable 
to the partially wrapped. Thus, both configurations can coexist over 
a wide range of droplet volume. At large membrane tensions (d), the 
partially wrapped configuration always has a lower energy. 
Nonetheless, the completely wrapped configuration is still 
metastable. (e-g) Typical morphologies with increasing droplet 
volume. 
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Fig. 4 Morphology diagram as a function of membrane tension 
and intrinsic contact angle with two distinct regimes: high and 
low tension regimes. The data points correspond to the full 
numerical solutions: circles indicate the occurance of budding 
transitions and crosses the absence of budding transitions. We also 
compare the numerical results with the simple prediction |T0| ), 
which is based on the spherical cap approximation. 



figurations are metastable over a wide range of droplet vol- 
umes, with the former having a lower energy than the latter 
within the computed volume range. It is, however, conceiv- 
able that if the completely wrapped configuration will become 
the state of lowest energy at an even larger droplet volume. In 
such a case, we compute the energy gradient dE/dV for the 
two morphologies. If this gradient is larger for the completely 
wrapped configuration in the limit of large droplet volumes, 
then we conclude that the system will not exhibit a budding 
transition. 



the completely and partially wrapped configurations represent 
the state of lowest energy for large and small volumes, respec- 
tively; and (II) a high tension regime, in which the partially 
wrapped configuration is the state of lowest energy for all vol- 
umes. The wetting properties of the liquids on the membrane 
play an important role. For smaller contact angles 0in, the 
liquid droplet wants to increase its contact area with the mem- 
brane, which favors the budded state. Therefore for smaller 
0in, the budding transition occurs over a larger range of mem- 
brane tension. 

We also note that, close to the boundary between regimes 
(I) and (II), it can be rather difficult to decide numerically if 
a budding transition is present. A typical example is provided 
by Fig. |3jd). Both the partially and completely wrapped con- 



5 Spherical cap approximation. 

We shall now argue that a simple, spherical cap approxima- 
tion as in Ref. may be used to estimate the boundary be- 
tween the two tension regimes shown in Fig. |4] Such an es- 
timate is rather useful since the full numerical solutions are 
time consuming to calculate. The spherical cap approxima- 
tion is based on the assumption that the contributions from the 
interfacial terms dominate the bending terms when the vol- 
ume of the droplet is large. This reflects the fact that the two 
energetic contributions scale differently with size. When this 
approximation is valid, the global shape of the system can be 
simply determined by considering the force balance of the in- 
terfacial tension with the membrane tensions. This is anal- 
ogous to Neumann's triangle for the force balance between 
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Fig. 5 Comparison of shapes as obtained from the full numerical 
solutions and from the spherical cap approximations for = 45°. 
The latter approximation breaks down at low membrane tension, but 
is accurate for high tension. The definitions of the tensions and 
effective angles are also shown. 



capillary surfaces'^. The resulting effective angles 6a, Oj^ and 
Oy, defined in Fig. [5j are given by 



cos 6i 



Y? - 



-I}.- 
ij 



(8) 



ij^ik 



with ij^k = a,j3,7 and / 7^ j k, compare also to Ref.'^'^^'. 
The sum of the effective angles is equal to 2;r. In contrast to 
capillary surfaces, for which Neumann's triangle applies, the 
membrane tensions, as well as the effective contact angles, do 
not represent material parameters; in contrast to the intrinsic 
contact angle as we have shown in^^. However, the three ten- 
sions are related via the relation for the intrinsic contact angle, 
Eq. ^. 

In Fig. [5] we compare some shapes as obtained from the 
spherical cap approximations with the full numerical solu- 
tions. At high membrane tensions, the spherical cap model 
predicts the partially wrapped configuration well, see Fig. 
[5ja). At lower membrane tensions, close to the boundary in 
the morphology diagram, the spherical cap approximation be- 
comes less accurate, compare Fig. [5jb), and it may give a so- 
lution corresponding to a partially wrapped configuration even 
when such a solution no longer exists for the original differ- 
ential equations described in section 2, see Fig. |5jc). Finally, 
at very low tensions, the spherical cap approximation does not 
have a solution anymore, as in Fig. |5jd). Using the force bal- 
ance at the contact line, the solution is lost when 



which implies 



^/3r ^ l+cosOir 

2 



(9) 



(10) 




Fig. 6 Evolution of membrane shape with increasing droplet volume 
for spontaneous curvature m = 5 ^Jt^^^^Jk, vanishing membrane 

tension E^y, and intrinsic contact angle Oi^ = 45°. The spontaneous 
curvature leads to nonspherical, tube-like membrane shapes. 



This restriction does not exist in the full numerical solutions 
because the bending terms can compensate the force imbal- 
ance between the interfacial and membrane tensions. This 
importance of bending terms in the limit of infinite droplet 
volume implies a small contact line radius, which is consis- 
tent with a droplet in the completely wrapped configuration. 
Therefore, the inequality ([TO]) represents a criterion for the 
possible existence of a completely wrapped state, and thus of a 
budding transition. The limiting case of relation ([TO]) as given 
by JLj^y = (1 + cos 0in)/2 corresponds to the dashed line in the 
morphology diagram. Fig. [4] Inspection of Fig. [4] shows that 
this line as obtained from the spherical cap approximation pro- 
vides a rather good estimate for the boundary between the low 
and high tension regimes as determined by the full numerical 
solution. Compared to the full numerical solutions, the spher- 
ical cap approximation always predicts that the budding tran- 
sition is lost at a lower membrane tension for a given intrinsic 
contact angle. 

In this study, we have ignored thermally excited undula- 
tions, which can be justified as follows. For a membrane ten- 
sion of the order of I^^p , the longest wavelength of these un- 
dulations is about (/c/Zof^)^/^. Using the typical experimen- 
tal valueslSl, ^ 10~^^J/m^ and K ^ 10~^ J, one obtains 
(/c/Lofjg)^/^ 2^ 100 nm, which is well below optical resolution. 
Furthermore, the area stored in these undulations is of the or- 



derof ll^ln 



times the membrane area 



10 



where /me 

is the small wavelength cut-off of the order of the membrane 
thickness. For /c 2^ 10 — 20 ksT at room temperature, the area 
fraction stored in thermally excited undulations is less than 
1-2 percents. 

6 Spontaneous curvature. 

So far, we have only considered the case where the membrane 
spontaneous curvature has a negligible effect. In this section, 
we shall briefly discuss the influence of the sponteneous cur- 
vature on the budding transition shown in Fig. [T] For small 
membrane tensions, this transition increases the bending en- 
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ergy of the membrane by about Stik, and decreases the interfa- 
cial energy by about ILg^^Va^ . The budding transition occurs 
when these two energies are comparable, i.e. for droplet size 
Va^ ~ (/c/Zofjg)^/^ 100 nm. The spontaneous curvature m 
should affect the budding transition when 1/m < {k/1.q^^Y^^ . 

This expectation is confirmed by our preliminary calcu- 
lations, see Fig. [5] Indeed, for m > (Zo^jg/K*)^/^, the re- 
sulting bud no longer has a spherical shape, but a tube-like 
structure. In fact, we find a variety of solutions to the dif- 
ferential equations presented in section 2. Detailed analysis 
on the global minimum configuration for a given liquid vol- 
ume is outside the scope of this paper and will be presented 
elsewhere. Nonetheless, it is interesting to note that we find 
tube-like membrane structures even in the absence of exter- 
nal pulling forces, as in Ref.^. The tubes are initiated by the 
minimization of the interfacial area Af^^g of the aqueous phases 
and their narrow structures are maintained by the spontaneous 
curvature of the membrane. 

7 Conclusions. 

In summary, partial wetting of membranes by liquid droplets 
leads to a new budding mechanism which is governed by the 
competition between bending and surface energies. Budding 
can be induced by increasing the volume of the liquid droplet, 
decreasing the membrane tension, and/or lowering the intrin- 
sic contact angle of the droplet at the membrane. Upon per- 
forming one of these changes, the membrane becomes more 
strongly bent until it undergoes a discontinuous transition to a 
completely wrapped configuration. 

Our theoretical results are accessible to experimental stud- 
ies. One example is provided by micropipette experiments, in 
which the budding transition may be induced by varying the 
membrane tension. Alternatively, by deflating the vesicles, a 
number of relevant parameters such as droplets volume, mem- 
brane tension and intrinsic contact angle can be altered at the 
same time to induce budding. 

In the experiments performed so far on these systems, 
the sizes of the two liquid compartments were comparable, 
whereas in our calculations, one phase is much larger than the 
other. It will be of great interest to extend our study and in- 
clude finite size effects arising from the finite volume of the j3 
phase. We expect that global constraints act to suppress bud- 
ding and only small buds will be possible. In this sense, the 
droplet-induced budding mechanism considered in this paper 
represents an important limiting behaviour. Other future ar- 
eas of interest include a detailed analysis of the effects arising 
from sponatenous curvature, and the influence of the wetting 
properties on the neck structure of the buds. 
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